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We present a general framework to study relativistic compound systems in a Hamiltonian formal- 
ism. This formalism is based on the explicitly covariant formulation of light-front dynamics, with 
a decomposition of the state vector in Fock components. In order to be able to make definite pre- 
dictions order by order in the truncation of the Fock expansion, we use an appropriate Fock sector 
dependent renormalization scheme. Our covariant formulation enables us to have a strict control 
of any violation of rotational invariance due to the choice of a given orientation of the light-front 
plane. This is mandatory in order to define the necessary renormalization conditions, and hence to 
be able to calculate physical observables. We emphasize the role played by antiparticle degrees of 
freedom in order to control, order by order in the Fock expansion, the scale invariance of physical 
observables. This nonperturbative framework demands also to use an appropriate regularization 
scheme. We show why the recently proposed Taylor-Lagrange regularization scheme is a very 
adequate scheme since it can be implemented very naturally, and from a systematic point of view, 
in light-front dynamics. As a direct application of this general framework, we settle the basis 
for a new formulation of chiral effective field theory in the baryonic sector based on light-front 
dynamics and a Fock decomposition of the state vector. 



X 



Light Cone 2010: Relativistic Hadronic and Particle Physics 
June 14-18, 2010 
Valencia, Spain 



* Speaker. 



© Copyright owned by the author{s) under the terms of the Creative Coinmons Attribution-NonCoinmercial-ShureAHke Licence. 



http: //pos. sissa. it/ 



Field theory on the light-front 



Jean-Frangois Mathiot 



The understanding of the structure of relativistic compound systems in nuclear and particle 
physics is the subject of intense theoretical studies since the discovery of the nucleon structure in 
the 60's. In both domains, the understanding of the properties of bound state systems from the 
original Lagrangian of quantum chromodynamics (QCD) or from an effective chiral Lagrangian, 
demands to develop a nonperturbative framework. The interest of light-front dynamics (LFD) to 
investigate relativistic bound states has been advocated for a long time. However, while its use 
in few-body systems has developed rapidly, its application to field theory beyond a perturbative 
framework is only a recent accomplishment. This is due to the fact that any practical calculation 
relies on the truncation of the Fock expansion of the state vector. 

We shall show in this study that we have now all the necessary theoretical tools in order to 
develop on a large scale a nonperturbative framework to calculate the structure of relativistic com- 
pound systems. Our formalism should enable us to control, order by order in the Fock expansion, 
the scale invariance of physical observables, and hence to make reliable predictions for these ob- 
servables. 

1. Covariant formulation of light-front dynamics 

LFD is one of the three forms of dynamics proposed by Dirac in 1949 [1] to describe physical 
systems, bound states as well as scattering states. In this form of dynamics, the state vector de- 
scribing the system is defined on the surface in four-dimensional space-time given by = ? + z/c. 
According to Dirac's classification, the ten generators of the Poincare group, given by space-time 
translations (four generators), space rotations (three generators), and Lorentz boosts (three genera- 
tors), can be separated into kinematical and dynamical operators. The kinematical operators leave 
the light-front plane invariant and are independent of dynamics, i.e. of the interaction Hamiltonian 
of the system, while the dynamical ones change the light- front position and depend therefore on the 
interaction. Among the kinematical operators, one finds, in LFD, the boost along the z axis. This 
property is of particular interest when one calculates electromagnetic observables at high momen- 
tum transfer. However, the spatial rotations in the xz and yz planes become dynamical, in contrast 
to the case of equal-time dynamics. This is a direct consequence of the violation of rotational in- 
variance caused by the non-invariant definition of the light-front plane orientation. This violation 
should be kept under control. 

While rotational invariance should be recovered automatically in any exact calculation, this is 
not a priori the case if the Fock expansion is truncated. In order to control the violation of rotational 
symmetry, we shall use the covariant formulation of LFD (CLFD) [2, 3], which provides a simple, 
practical, and very powerful tool in order to describe physical systems. In this formulation, the 
state vector is defined on the plane characterized by the invariant equation (O x = a, where (O is 
an arbitrary light-like (ft)^ = 0) four-vector, and a is the Ught-front time. The standard LFD is 
recovered by considering the particular choice CD = (1,0,0, —1). The covariance of our approach 
relies on the invariance of the light-front plane equation under any Lorentz transformation of both 
(O and X. This implies in particular that ft) cannot be kept the same in any reference frame, as it is 
the case in the standard formulation of LFD. 

There is of course equivalence between the standard and covariant forms of LFD in any exact 
calculation. Calculated physical observables must coincide in both approaches, though their deriva- 
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tion in CLFD in most cases is much simpler and more transparent. Indeed, the relation between 
CLFD and standard LFD reminds that between the Feynman graph technique and old-fashioned 
perturbation theory. In approximate calculations, CLFD has a definite advantage in the sense that 
it enables a direct handle on the contributions which violate rotational invariance. These ones de- 
pend explicitly on the orientation of the light-front surface (i.e. on (o) and can thus be separated 
covariantly from true physical contributions. 

The transformation properties of the light-front plane are governed by the four-dimensional 
momentum and angular momentum operators given by 

= po + p^"' with P'^' = (0^1 H'"' {x)d{(o-x-a)d^x , (1.1) 

f^iv = J^y +J^llv with = j H'"' {x){x^(Ov -xy(0^)5{(0-x - o)d^x . (1.2) 

In these equations, refers to the interaction Hamiltonian, while (-/JSy) ^^e standard free 
four-momentum (four-angular momentum) operator. According to the general properties of LFD, 
we decompose the state vector <^{p) of a physical system of momentum p in Fock sectors, with 

oo p 

^ip)=J^ dDn<l>„{ki,...,kn;p)5'^{ki + ...+kn-p-(0Tr,)\n) , (1.3) 

«=!"' 

where \n) is the state containing n free particles with four- momenta k\,...,kn and (^,,'s are rela- 
tivistic n-body wave functions, the so-called Fock components. The phase space volume is de- 
noted by dDn. Note the particular overall momentum conservation law given by the 5-f unction. 
It follows from the general transformation properties of the light-front plane co x = G under the 
four-dimensional translations governed by in (1.1). The quantity T„ is a measure of how far the 
«-body system is off the energy shell. It is completely determined by this conservation law and 
the on-mass-shell condition for each individual particle momentum. It is convenient to introduce, 
instead of the wave functions 0„, the vertex functions r„ defined by r„ = {s„ — M^)(^„, where M is 
the physical mass of the bound state and 5,, = + ^2 + • • O^- The system of coupled equations for 
the Fock components of the state vector can be obtained from the general eigenstate equation [4] 

P'0(p)=m20(p). (1.4) 

The state vector is finally normalized to 1. 

As follows from the transformation properties of the state vector under the four-dimensional 
angular momentum (1.2) (the so-called angular condition [3]), the spin structure of the wave func- 
tions <p„ is very simple, since its construction does not require the knowledge of dynamics. It 
should incorporate however G)-dependent components. It is convenient to decompose each wave 
function 0„ into invariant amplitudes constructed from all available particle four-momenta (includ- 
ing the four- vector w!) and spin structures (matrices, bispinors, etc.). In the Yukawa model for 
instance (one constituent fermion coupled to scalar bosons), we have for the one- and two-body 
components [4] 



Ti = aiu{ki)u{p), r2 = u{ki 



bi+bj 

cop 



u{p), (1.5) 



since no other independent spin structures can be constructed. Here ai, bi, and Z72 ai^e scalar 
functions determined by dynamics, and ^1 is the momentum of the constituent fermion. 
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2. The Fock sector dependent renormalization scheme 

The application of LFD to bound state systems in field theory necessitates to define an appro- 
priate, nonperturbative, renormalization scheme. Indeed, the full cancellation of field-theoretical 
divergences which appear in a given Fock sector requires to take into account contributions from 
other sectors. If even a part of the latter is beyond our approximation, some divergences may leave 
uncancelled. This is illustrated in Fig. 1 for the calculation of the fermion propagator in second 
order of perturbation theory: the self-energy contribution and the fermion mass counterterm in- 
volves two different Fock sectors. This means that, as a necessary condition for the cancellation of 
divergences, any mass counterterm should be associated with the number of particles present (or 
"in flight") in a given Fock sector. In other words, all mass counterterms must depend on the Fock 
sector under consideration, as advocated first in Ref. [5]. This is also true for the renormalization 
of the bare coupling constant. 



Figure 1: Renormalization of the fermion propagator in the second order of perturbation theory. The fermion 
mass counterterm is denoted by 5m. 

The presence of uncancelled divergences reflects itself in possible dependence of (approx- 
imately) calculated observables on the regularization parameters, i.e. approximate physical ob- 
servables are not anymore scale invariant. This prevents to make any physical predictions if we 
cannot control the renormalization procedure in one way or another. To do that, we have developed 
an appropriate renormalization procedure — the so-called Fock sector dependent renormalization 
(FSDR) scheme — in order to keep the cancellation of field-theoretical divergences under per- 
manent control. This scheme amounts to consider a complete set of mass counterterms and bare 
coupling constants, denoted by 5m^"^ and ^q"' respectively, where n refers to the number of bosons 
in flight [4]. The mass counterterm is fixed from the solution of the eigenstate equation (1.4), by 
demanding that the mass of the constituent fermion is identical to the mass of the physical one. The 
bare coupling constant is fixed from the general, on energy shell, renormalization condition [6] 



where Zf and Zi, are the field strength renormalization factors of the constituent fermion and boson 
respectively. When the Fock space is truncated, one should also make sure that all contributions 
remain within the Fock space. 

The condition (2.1) has two important consequences. The first one is that the two-body vertex 
function at S2 = should be independent of the four- vector ft) which determines the orientation 
of the light- front plane. With the spin decomposition (1.5), this implies that the component b2 at 
S2 = should be identically zero. While this property is automatically verified in the case of the 
two-body Fock space truncation - if using a regularization scheme which does not violate rotational 
invariance - this is not guaranteed for calculations within higher order truncations. Indeed, nothing 
prevents r2 to be ftj-dependent, since it is an off-shell object, but this dependence must completely 
disappear on the energy shell, i.e. for S2 = M^. The truncation of the Fock space may cause 




+ 



Srn 




(2.1) 
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some ft)-dependence of T2 even on the energy shell, which makes the general renormalization 
condition (2.1) ambiguous. If so, one has to insert new counterterms into the light-front interaction 
Hamiltonian, which explicitly depend on (O and cancel the spurious ft)-dependence of r2(52 =M^). 
Note that the explicit covariance of CLFD allows to separate the terms which depend on the light- 
front plane orientation (i.e. on ft)) from physical contributions, and establish the structure of these 
counterterms. This is not possible in ordinary LFD. 

The second consequence is that r2{s2 = M^) should be independent of x. This is a non- 
trivial requirement since r2 does depend in general on two invariant kinematical variables which 
are usually chosen as the longitudinal momentum fraction, x, of the constituent fermion, and the 
square of its transverse momentum, k^. Hence, if we fix from the on energy shell condition 
S2 = M^, r2 should be independent of x. Again, this property is verified in the two-body Fock space 
truncation, since our equations are equivalent in this approximation to perturbation theory of order 
It is not guaranteed for higher order calculations. In practice, we shall fix r2{s2 = M^) at some 
preset value x* and verify that calculated physical observables are not sensitive to the choice of x*. 

We show in Figs. 2 the anomalous magnetic moment in the Yukawa model [6], in the three- 
body Fock state truncation, for two typical values of the boson-fermion coupling constant. These 
calculations do not involve antiparticle degrees of freedom. At relatively small values of the cou- 
pling constant, they show a rather nice scale invariance, while the calculation at higher values of 
the coupling constant shows the onset of scale dependence. 
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Figure 2: The anomalous magnetic moment in the Yukawa model as a function of the regularization mass /ii 
(Pauli-Villars boson mass), for two different values of the coupling constant, a = 0.2 (left plot) and 0.5 (right 
plot). The dashed and dotted lines are, respectively, the two- and three-body contributions, while the solid 
line is the total result. The value of the anomalous magnetic moment calculated in the N — 2 approximation 
is shown by the thin line on the right axis. 

We identify this (small) residual scale dependence with the lack of contributions from antipar- 
ticle degrees of freedom which should be incorporated explicitly in LFD. To have a better handle on 
these contributions, it is instructive to calculate ^1.2(^2 = M^) in perturbation theory. This is done 
by calculating the amplitudes of the diagrams shown in Fig. 3. The contribution of the left diagram 
in Fig. 3 is automatically taken into account by the solution of the eigenstate equation in the N = 3 
Fock state truncation, while the right plot corresponds to the contribution of fermion-antifermion 
contributions in leading order. If one calculates the sum of both contributions [7], one finds 

bP'"''{s2=M^) = cte , (2.2) 
bP"''{s2=M^) = 0. (2.3) 
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This is a first indication that the expected properties of the on-shell functions b\^2 are indeed recov- 
ered when antifermion degrees of freedom are involved. 




Figure 3: Part of the three-body Fock sector contributions to the two-body vertex. The thick solid line 
represents an antifermion. 

3. The Taylor-Lagrange regularization scheme 

Before choosing an appropriate regularization scheme, it is interesting to come back to the 
very origin of divergences of loop integrals. These divergences can be traced back to the violation 
of causality, originating from ill defined products of distributions at the same point [8]. The correct 
mathematical treatment, known since a long time [9], is to consider covariant fields as operator 
valued distributions (OPVD), these distributions being applied on test functions with well-defined 
properties. These considerations lead to the Taylor-Lagrange renormalization scheme (TLRS) [10]. 

If we denote by / the Fourier transform of the test function, a scalar- field <p(x), for instance, 
will thus write 

with £p = -\-m^. From this decomposition, it is apparent that test functions should be attached 
to each fermion and boson fields, while the contraction of two fields (propagator) should be pro- 
portional to In LFD, / depends only on p^. These test functions should have three important 
properties: 

/) the physical field (p{x) in (3. 1) should be independent of the choice of the test function. This 
may be achieved if / is chosen among the partitions of unity (PU), i.e. if / is build up from a family 
of functions j3, (X) with 

L 

£A(^) = 1 for X^[a,b\ . (3.2) 
/=i 

It is a function of finite support which is 1 everywhere except at the boundaries. This choice is also 
necessary in order to satisfy Poincare invariance [10]. 

//) In order to be able to treat in a generic way singular distributions of any type, the test 
function is chosen as a super regular test function (SRTF), i.e. a function which vanishes, as well 
as all its derivatives, at the ultraviolet (UV) and the infrared (IR) boundaries. 

///) The boundary conditions of the test function - which is assumed for simplicity to depend 
on a one dimensional variable X - should embody a scale invariance inherent, in the UV domain 
for instance, to the limit X — )• oo since in this limit rj^X also goes to oo, where T]^ is an arbitrary 
dimensionless scale. This can be done by considering a running boundary condition for the test 
function, i.e. a boundary condition which depends on the variable X according to 

fix > H{X)) = for H{X) = 71^Xg{X) . (3.3) 
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This condition defines a maximal value, X^ax, with f{Xmax) = 0. A typical example of the 
fuction g{X) is given in Ref. [10]. This running condition is equivalent to having an ultra-soft 
cut-off, i.e. an infinitesimal drop-off of the test function in the asymptotic limit, the rate of drop- 
off being governed by the arbitrary scale rj^. A similar scale invariance is also present in the IR 
domain. An example of test function constructed from two elementary functions, with a running 
condition, is shown on Fig. 4 [10] 




Figure 4: Construction of a partition of unity with running support. The left curve shows the IR domain, 
while the right curve shows the UV domain. 

With these properties, the TLRS can be summarized as follows, in the UV domain for instance. 
Starting from a general amplitude £/ written for simplicity in a one dimensional space as 

poo 

£/= dX T{X) f{X) , (3.4) 
Jo 

where r(X) is a singular distribution in the UV domain, we apply the following general Lagrange 
formula to f{X), after separating out an intrinsic scale a from the (running) dynamical variable X 



(3.5) 



This Lagrange formula is valid for any order k, with k> Q, since / is chosen as a SRTF and it is 
therefore equal to its Taylor remainder for any k. After integration by part in (3.4), and using (3.5), 
we can express the amplitude jz/ as 

/>oo 

dX T{X) f{X) , (3.6) 





where r(X) is the so-called extension of the singular distribution T{X). In the limit / —)• 1 obtained 
by letting — ;> 1^, it is given by [10] 

The value of k in (3.7) corresponds to the order of singularity of the original distribution T{X) 
[10]. The extension of T{X) is no longer singular due to the derivatives in (3.7), so that we can 
safely perform the limit / — )• 1 in (3.6). This amplitude is well defined but depends on the arbitrary 
dimensionless scale T]^. This scale is the only remnant of the presence of the test function. Note 
that we do not need to know the explicit form of the test function in the derivation of the extended 
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distribution f{X). We only rely on its mathematical properties and on the running construction 
of the boundary conditions. The extension of singular distributions in the IR domain can be done 
similarly. 

The use of TLRS for the calculation of the state vector of compound systems within CLFD is 
very natural. Since each vertex function r„ is attached to one fermion and n — \ boson lines, it will 
be multiphed at least by /(k^) . . ./(k^). We can thus redefine r„ to include implicitly these test 
functions: 

r„(ki . . .k„) ^ r„(ki . . .k„) = r„(ki . . .k„)/(k?) . . ./(k^) . (3.8) 

Since the f's are SRTF, this implies that any r„ is also a super regular function with respect to 
all momenta. We can thus apply the Lagrange formula to all loop calculations, and derive the 
extension of all singular distributions along the lines detailed in the previous section [10]. 

It is instructive to calculate the extension of the singular distribution T{X) = \/{X + a). Using 
(3.5) and (3.6) with k = 0, one finally aiiives to the following extension [11] 

fix) = — ^-^ . (3.9) 

This form is reminiscent of a Pauli-Villars subtraction, with however two important differences. 
This subtraction involves the arbitrary scale rf'. We thus do not have to perform any infinite mass 
limit. Moreover, this extension results from the application of TLRS to physical contributions, and 
we do not have to introduce any ad-hoc Pauli-Villars particles. 



4. Light front chiral effective field theory 

It is now well recognized that the properties of QCD at low energies can be understood in terms 
of an effective Lagrangian based on the (approximate) chiral symmetry of strong interactions, and 
its spontaneous breakdown. The physical properties of pions within chiral perturbation theory 
(CPT) are now well reproduced in a meaningfull power expansion. This originates from the fact 
that the pion has zero mass in the chiral limit. This is however not the case in the nucleon sector, 
and all momentum scales are a-priori involved in the calculation of nucleon properties. 

Since the pion mass is zero in the chiral limit, any calculation of nN systems demands a 
relativistic framework to get, for instance, the right analytical properties of the physical amplitudes. 
The calculation of compound systems, like a physical nucleon composed of a bare nucleon coupled 
to many pions, relies also on a nonperturbative eigenstate equation. While the mass of the system 
can be determined in leading order from the iteration of the nN self-energy calculated in the first 
order of perturbation theory, as indicated in Fig. 5(a), this is in general not possible, in particular, 
for TiN iiTcducible contributions, as shown in Fig. 5(b). 

While there is not much freedom, thanks to chiral symmetry, for the construction of the effec- 
tive Lagrangian in CPT in terms of the pion field n — or more precisely in terms of the U field 
defined hy U = e^'^-^lf", where fj^ is the pion decay constant and T are the Pauli matrices, — one 
should settle an appropriate approximation scheme in order to calculate nucleon properties. Up to 
now, two main strategies have been adopted. The first one is to force the bare (and hence the phys- 
ical) nucleon mass to be infinite, in heavy baiyon chiral perturbation theory [12]. In this case, by 
construction, an expansion in characteristic momenta can be developed. The second one is to use 
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1 \—k \ 1^-4 \^ 

(a) (b) 

Figure 5: Iteration of the self-energy contribution in the first order of perturbation theory (a); irreducible 
contribution to the bound state equation (b). Dashed lines represent pions. 

a specific regularization scheme [13] in order to separate contributions which exhibit a meaningful 
power expansion, and hide the other parts in appropriate counterterms. In both cases however, the 
explicit calculation of baryon properties relies on an extra approximation in the sense that physical 
amplitudes are further calculated by expanding the effective Lagrangian, denoted by ^cpt, in a 
finite number of pion fields in perturbation theory. 

The general framework we have developed above is particularly suited to deal with these 
requirements in the nucleon sector. This leads to the formulation of light-front chiral effective field 
theory (LF;^EFT) [14] with a specific effective Lagrangian ^eff- The decomposition of the state 
vector in a finite number of Fock components implies to consider an effective Lagrangian which 
enables all possible elementary couplings between the pion and nucleon fields compatible with 
the Fock space truncation. This is indeed easy to achieve in CPT since each derivative of the U 
field involves one derivative of the pion field. In the chiral limit, the chiral effective Lagrangian of 
order p, ^Qpj, involves p derivatives and therefore at least p degrees of the pion field. In order to 
calculate the state vector in the A/^-body approximation, with one fermion and (A^ — 1) pions, one 
has therefore to include contributions up to 2 (A'^ — 1) pion fields in the effective Lagrangianv (for 
pions in the initial and final states). We thus should calculate the state vector in the A'^-body Fock 
state truncation with an effective Lagrangian denoted by -^^ff with 

^^fj = ^i;?^-'\ (4.1) 

While the effective Lagrangian in LF;^EFT can be mapped out to the CPT Lagrangian of order 
p, the calculation of the state vector does not rely on any momentum decomposition. It relies only 
on an expansion in the number of pions in flight at a given light-front time. In other words, it 
relies on an expansion in the fluctuation time, if, of such contributions. From general arguments, 
the more particles we have at a given light-front time, the smaller the fluctuation time is. At low 
energies, when all processes have characteristic interaction times larger than T/, this expansion 
should be meaningful. 

Moreover, the use of TLRS enables to perform systematic nonperturbative calculations in the 
general framework developped in Sees. 1 and 2, contrarily to the regularization schemes used up 
to now. This finite regularization scheme is also of particular interest since it does not involve any 
large momentum or energy scale which may break chiral symetry. 

It is interesting to illustrate the general features of LF;^EFT calculations. At order A'^ = 2, we 
already have to deal with irreducible contributions, as shown in Fig. 5(b). It leads to a non-trivial 
renormalization of the coupling constant. The calculation at order N = 3 incorporates explicitly 
contributions coming from titi interactions (a and p resonances), as well as low energy TiN res- 
onances, like the A or Roper resonances. We can generate therefore all these resonances in the 
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intermediate state without the need to include them explicitly, provided the effective Lagrangian 
has the right dynamics to generate these resonances. 

5. Perspectives 

The general framework we presented above to study the structure of relativistic compound 
systems in light-front dynamics in a nonperturbative way relies on three main advances: i) the con- 
struction of a covariant formulation of light-front dynamics [2, 3] in order to control any violation 
of rotational invariance; ii) the development of an appropriate renormalization scheme — the so- 
called Fock sector dependent renormalization scheme — to deal with the truncation of the Fock 
expansion [4] ; Hi) the use of an appropriate regularization scheme — the so-called Taylor-Lagrange 
regularization scheme — very well adapted to systematic calculations in light-front dynamics [11]. 

These advances should enable us to have a predictive framework order by order in the Fock 
expansion. We shall complete in the future this description by considering physical systems involv- 
ing spontaneous symmetry breaking. It is known that these systems can be described in light-front 
dynamics by the consideration of zero-mode contributions, in the theory in 1 + 1 dimension 
for instance [ 15] . Their full calculation in 3 + 1 dimensions within the general framework presented 
above remains to be done. 
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